Noncrossing Approximation for the Anisotropic Kondo Model: Charge Fluctuations 

in a Quantum Box 
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The noncrossing approximation (NCA) is generalized to the multi-channel Kondo-spin Hamilto- 
nian with arbitrary anisotropic exchange couplings and an external magnetic field, and applied — 
in the framework of Matveev's mapping — to the charge fluctuations in a single-electron box at 
the Coulomb blockade. The temperature dependences of the charge step and the capacitance are 
calculated for a narrow point contact. At low temperatures and close to the degeneracy point, the 
capacitance line shape exhibits an approximate scaling with U /VT, where U is the deviation in 
gate voltage from the degeneracy point. This scaling relation is proposed as a sharp experimental 
diagnostic for the non-Fermi-liquid physics of the system at low temperatures. Both the reliability 
and shortcomings of the Kondo NCA are discussed in detail. Through comparison with poor-man's 
scaling, we are able to pinpoint the omission of particle-particle processes as the origin of the NCA 
flaws. An extended diagrammatic scheme is devised to amend the NCA flaws. 



PACS numbers: 72.15.Qm, 73.23.Hk, 75.20.Hr 
I. INTRODUCTION 

The rapid development of nanofabrication techniques 
has opened new possibilities for the study of the Kondo 
effect in mesoscopic systems. Using ultra small quantum 
dots it is now possible to measure Kondo-assistedlunncl- 
ing through a single tunable magnetic impurityuU Using 
scanning tunneling spectroscopy, one can directly probe 
the local electronic structure around an isolated magnetic 
adatom on a metallic surfacefrtl It is even becoming ex- 
perimentally feasible to address the competition between 
intersite magnetic locking and Kondo-type correlations 
in the two-impurity cluster .EfE3 Perhaps most intriguing, 
though, arc tJaft-eossible realizations of the two-channel 
Kondo effect, Ou3 which is a prototype for a non-Eermi- 
liquid ground state in correlated electron systems .L3 

There are several leading candidates for the observa- 
tion of the two-channel Kondo effect in mesoscopic sys- 
tems. Most notabhz, .zero-bias anomalies seen in metal- 
lic point contacts,cj~t3 which feature square-root tem- 
perature and voltage dependence of the differential con- 
ductance; scaling of the differential conductance with 
eV/ksT (eV being the applied bias); and the absence of a 
Zeeman splitting under an applied magnetic field — all of 
which are in line with the notion of two-channel Kondo 
scattering off nonmagnetic two-level tunneling systems 
(TLS) Mill Two-channel Kondo scattering off TLS was 
also recently invokedta to explain the shape and scal- 
ing behavior of the out-of-equilibrium energy distribution 
function of quasiparticles in copper and gold wires JIB In 
a separate recent experiment on the changing of a semi- 
conductor quantum dot, Berman et alEU observed ca- 
pacitance line shapes in accordance with Matveev's pre- 
diction for the two-channel Kaadq screening of quantum 
charge fluctuations on the dotHEl The two nearly degen- 
erate levels in this picture correspond to the two available 
charge configurations at the Coulomb blockade. 



Common to both the TLS and Coulomb-blockade sce- 
narios for the two-channel Kondo effect is a strong 
anisotropy in the effective spin-exchange interaction. 
Specifically, in TLS one generically hasEJ J z >• J x > 
J y = (J, being the ith component of the effective 
spin-exchange interaction), while J z — and—J^ = 
J y 7^ in the case of the Coulomb blockaded That 
the intermediate-coupling fixed point of the two-channel 
Kondo model is stable against exchange anisotrony-is. 
well known from numerical penormalization-groupElL-j 
and conformal field theoryniiEa studies. However, there 
is a clear lack in quantitative treatments that cover all 
temperature and field regimes of the model, as is re- 
quired for each of the two systems mentioned above. 
This is especially true of the nonequilibrium scattering off 
TLS, where theoretical efforts have focused-on the related 
SU(2)xSU(2) two-channel Anderson modelJlJ rather than 
the actual anisotropic Kondo-spin Hamiltonian appropri- 
ate for this case. We also note a growing interest in quan- 
titative-theories of the single-channel anisotropic Kondo 
model,E3 particularly in.c<j)ftiiection with the Ohmic dissi- 
pative two-state systep. — ciil and the heavy fermion com- 
pound Ceo.8Lao.2Al3.EtE 

In—this paper, we extend the noncrossing approxima- 
tion^ (NCA) to the Kondo -spin Hamiltonian with arbi- 
trary anisotropic spin-exchange couplings, and apply it 
to the charge fluctuations in a single-electron box at the 
Coulomb blockade. The NCA was extensively used in the 
80's to study dilute magnetic alloys, especially Ce- and 
Yb-based alloys with large orbital degeneracy. Its useful- 
ness for treating the multi-channel Andersan model was 
later emphasized by Cox and Ruckenstein,E3 who noticed 
that—the NCA pathology hampering the single-channel 
caseEl actually corresponds to the exact non- Fermi-liquid 
power laws and logarithms of the multi-channel Kondo 
effect. Compared to exact methods such as the Bcthe 
ansatz and conformal field theory, the NCA has the cru- 
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cial advantage that it is not restricted to idealized mod- 
els; it could be used to compute dynamical properties 
over extended temperature and parameter regimes; and 
it has a natural extension to nonequilibriumO'Lj Aside 
from the Anderson model, the NCA was also applied to 
the Coqblin-Schrieffcr Hamiltoniai£3 in the form of a self- 
consistent ladder approximation,!^ yet no direct formu- 
lation of this approach has been devised to date for the 
Kondo Hamiltonian. 

Generalizing the NCA to the multi-channel Kondo- 
spin Hamiltonian with arbitrary spin-exchange and 
potential-scattering couplings, a comprehensive analysis 
of this approach is presented. It is shown that the Kondo- 
NCA (KNCA) correctly describes the non-Fermi-liquid 
fixed point of the multi-channel Kondo effect, similar 
to the-MCA formulation of the multi-channel Anderson 
model.tS Furthermore— there is quantitative agreement 
with the Bethe ansatzEa for the temperature and field de- 
pendence of the magnetic susceptibility of the isotropic 
two-channel model. At the same time, the KNCA has 
several shortcomings, including a spurious ferromagnetic 
Kondo effect, and an inaccurate exponential dependence 
of the Kondo temperature on the inverse coupling con- 
stants. We are able to pinpoint these shortcomings of 
the NCA as due to the omission of particle-particle di- 
agrams, which also seems to be the reason for the well- 
known NCA failure to describe the Fermi-liquid ground 
state of the single-channel Kondo effect. As a first step 
towards amending this flaw, an extension of the KNCA is 
constructed, incorporating particle-particle and particle- 
hole diagrams on equal footing. 

Applying the KNCA to the charge fluctuations in a 
quantum box with a single-mode junction, we propose 
an approximate scaling relation for the capacitance line 
shape as a sharp experimental diagnostic for the low- 
temperature, non- Fermi- liquid regime of the two-channel 
Kondo effect. This approximate scaling relation should 
prove useful in analyzing future experiments, especially 
on metallic boxes which are more favorable to display a 
fully developed two-channel Kondo effect .e2I 

The remainder of the paper is organized as follows: 
Section || briefly reviews Matveev's mapping of the 
charge fluctuations in a single-electron box onto an 
anisotropic two-channel Kondo Hamiltonian. The KNCA 
is then presented and analyzed in Sec. III. Technical de- 



tails of the KNCA and an extended discussion of aspects 
of this approach not directly related to the Coulomb 
blockade are deferred to Appendixes |a| and [§]. The ap- 
plication of the KNCA to the charge fluctuations at the 
Coulomb blockade is reported in turn in Sec. |v|, followed 
by a discussion of the main results in Sec. |v|. Appendix ^ 
is devoted an extension of the KNCA that properly pre- 
serves particle-hole symmetry. 



II. MAPPING THE COULOMB BLOCKADE 
ONTO AN ANISOTROPIC KONDO PROBLEM 

The fundamental connection between the charge fluc- 
tuations in a single-electron box at the Coulomb blockade 
and the Kondo problem was formulated in a key paper by 
MatveevEil Below we briefly outline Matveev's mapping 
of the two problems, setting our notations in doing so. 

Consider a single-electron box, connected to a lead by 
a narrow point contact. For a sufficiently narrow point 
contact, there will be just a single mode weakly connect- 
ing the two sides of the constriction. Such a weak link 
can be described by the tunneling Hamiltonian 
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where c ka (c^ CT ) creates a conduction electron with spin 
projection a in the lead (box), and tk p are the corre- 
sponding tunneling matrix elements. To simplify our 
notations, we use the indices k and p throughout this 
section to label the single-electron levels in the lead and 
in the box, respectively. We further assume that the level 
spacing inside the quantum box is sufficiently small that 
a continuum- limit descriptionj-ean be used. As recently 
emphasized by Zarand et a/.,c3 it is difficult to realize 
this condition in present semiconducting devices, due to 
the exponentially small many-body energy scale (Kondo 
temperature) involved. We shall return discuss this as- 
sumption later on. 

In addition to the tunneling Hamiltonian of Eq. (Q), 
one has the terms describing the isolated lead and box. 
While the lead can be modeled by a simple noninteracting 
band with dispersion efc, 
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for the quantum box one has to consider also the charging 
energy of the box. Specifically, setting the filled Fermi sea 
as our reference state and measuring all single-particle 
energies relative to the Fermi level, the charge inside the 
box is given by 



which has the associated charging energy 

Q 2 _ 
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Here C'b is the capacitance of the quantum box, and <p is 
the electrostatic potential in the box measured relative to 
the lead. The latter is related to the applied gate voltage 
by some geometric factor. Thus, the full Hamiltonian of 
the system has the form H = Hl + Hb + Ht, where 
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FIG. 1. Zero-temperature charging of the isolated quan- 
tum box, in the absence of tunneling to the lead. For a given 
value of the electrostatic potential <f), the charging energy 
Ec = (eN) 2 /2Cb — e<f>N is parabolic in the number of ex- 
cess electrons in the box, TV = —Q/e. Each time <^>Cs/e is 
equal to half an integer, there are two degenerate charge con- 
figurations that minimize Ec- This results in abrupt jumps 
in N each time cj>CB/e crosses a half-integer value. 



describes the isolated single-electron box, which has the 
single-particle levels e p . 

For zero tunneling and zero temperature, the charge 
in the quantum box is determined by minimizing Ec 
with respect to Q. Given the discreteness of the charge, 
this results in abrupt changes in the ground state of the 
quantum box each time 4>Cb /e crosses the half-integer 
value N + 1/2, corresponding to the transition from 
Q = -eN to Q = -e(N + 1) (see Fig. |). Precisely 
at tpCs/e = N + 1/2 the two charge configurations are 
degenerate, giving rise to strong charge fluctuations in 
the box. 

The effect of weak tunneling is to smear the sharp 
charge steps at <f> = e(N + ^)/Cb, yet in a nontrivial 
manner to be explored below. To this end, let us fo- 
cus on the vicinity of a particular charge step separating 
Q = —eN from Q = —e(N +1). Labeling the deviation 
from the degeneracy point by U = (f> — e(N + \)/Cb, 
we concentrate on \U\ <C e/Cs and ksT <C c 2 /Cb, such 
that all charge configurations other than Q,— —eN and 
Q = —e(N + 1) are thermally inaccessible o One can re- 
move all higher energy charge configurations in this case 
by introducing the operators Pq and Pi, which project 
onto the Q = —eN and Q = —e(N + 1) subspaces, re- 
spectively. Omitting a [/-independent reference energy, 
the resulting low-energy Hamiltonian is written as 
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eff 
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\ ka pa / 

+ X {tkpPocl^CpaPi + tlpPiC^CkaPoj 
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To illustrate the connection between Eq. (|fj) and the 
two-channel Kondo Hamiltonian, it is useful to adopt the 
following 2x2 matrix representation of the two projection 
operators: 



Po 
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Pi 
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(7) 



In addition, we change our notation for the conduction 
electrons by introducing the isospin label a = ±, distin- 
guishing the quantum-box electrons, c s kaa= _, from the 



lead electrons, 



Neglecting all momentum de- 



pendence of the tunneling matrix elements, i.e., taking 
tkp — t, Eq. (^|) acquires the form 
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Finally, expressing the 2x2 matrices in terms of a spin-i 
isospin operator S, our effective Hamiltonian becomes 

n eff = J2 e kA aa c kaa - eUS z - eU{N + 1/2) 

kaa. 

+ l X { C kaa a aa' C k'a a 'S + + h.c} , (9) 
kk' acta' 

where = a x ± ia y are Pauli matrices. Here we 
have omitted for convenience the isospin index from the 
conduction-electron energies in the first term of Eq. (JTJ), 
which amounts to taking a single joint density of states, 
p(e) = J2k — £ fe)' f° r t ne quantum box and the lead. 
Since all quantities related to S in the problem depend on 
the product of the two density of states in the box and 
in the lead, p{e) is to be understood as the geometric 
average of the two density of states .l3 

Up to a constant shift in energy, Eq. (|TJ) is the Hamil- 
tonian of an anisotropic two-channel Kondo impurity in 
a local magnetic field. Specifically, J± = 2t plays the role 
of the transverse exchange coupling, while pBdjH = sU 
corresponds to an applied magnetic field. The longitu- 
dinal exchange coupling in this mapping is set equal to 
zero. Note that the exchange interaction in Eq. ([)]) acts 
on the isospin index, whereas the physical spin is reduced 
to a passive spectator. Hence the associated many-body 
screening is that of the charge fluctuations in the quan- 
tum box, not that of a physical spin. More specific, the 
charge in the box is related to the isospin magnetization 
through 
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while the capacitance of the junction, C(U, T) 
—d(Q) /SLLis equal to e 2 / (pb9j) times the isospin 
ceptibility,a 



sus- 



X{H,T) = fi B gj 



d(Sz) 
dH 



(11) 



As is well known, the low-temperature, low-field prop- 
erties of the two-channel Kondo Hamiltonian are gov- 
ernedJay an intermediate-coupling non- Fermi-liquid fixed 
point J13 which differs markedly from the infinite-coupling 
fixed point of the single-channel case. In particular, the 
zero-field susceptibility diverges with decreasing temper- 
ature according to x cx (1/Tr-) log(Tjf/T), which trans- 
lates in the present context to an infinite charging slope 
at U — and T = 0. T# in the above expression 
is the Kondo temperature, which marks the crossover 
from weak coupling. pEpr the anisotropic Hamiltonian of 
Eq. (|) it is given byfeil 



k B T K = (Dp t) exp 



A Po t 



(12) 



where D ~ £ 2 /Cb is the effective conduction-electron 
bandwidth, and po is the geometric average of the lead 
and the box density of states at the Fermi level. 

Equation ([l2|) is due to Matveev, who studied 
the Hamiltonian of Eq. (^|) within a perturbative 
renormalization-group approach.Eil Based on perturba- 
tion theory, this approach breaks down as soon as the 
renormalized exchange couplings are large, limiting its 
applicability to the weak-coupling regime. Our objective 
is to develop a quantitative description of all temper- 
ature and field regimes, including the low-temperature, 
non-Fermi-liquid regime. To this end, we first extend 
the noncrossing approximation to the Kondo Hamilto- 
nian with arbitrary anisotropic spin-exchange couplings. 



III. NONCROSSING APPROXIMATION 

As mentioned above, to date there is no direct formu- 
lation of the NCA for the Kondo-spin Hamiltonian, let 
alone its strongly anisotropic version of Eq. (^). In fact, 
previous applications of the NCA to the multi-channel 
Kondo effect have been restricted to the SU(iV)xSU(M) 
multi-channel Anderson model.Ej While the two-channel 
Anderson model has often been used to model uranium 
and cerium ions in appropriate crystalline-electric-field 
environments, c3 it does not directly apply to the prob- 
lem at hand, or to any realization of two-level systems 
in a metal. Furthermore, we note that the anisotropic 
Kondo Hamiltonian cannot be derived from the Anderson. 
Hamiltonian via the Schrieffer- Wolff transformation,!^ 
which illustrates the independence of the two models. 
Below we extend the formulation of the NCA to the 
multi-channel Kondo Hamiltonian with arbitrary spin- 
exchange and potential-scattering couplings. 
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FIG. 2. Diagrammatic representation of the noncrossing 
approximation for the pseudo-fermion Green function. Here 
thin (thick) dashed lines represent bare (fully dressed) / 
Green functions, while full lines stand for bare conduc- 
tion-electron propagators. In the Coqblin-Schrieffer limit, 
Jo — Jx = J y — Jz, these diagrams coincide with t]*« 
self-consistent ladder approximation of Maekawa et 
which can also be derived from the NCA equations for the 
Andofspn impurity model by taking the appropriate Kondo 
limitEfl 



As we are interested in devising a general approach, 
let us consider in this section the fully anisotropic M- 
channcl Kondo Hamiltonian given by 

M 

K = YsYs X] £kcl na c knrT - p B gjHS z 
k n=i o=T 4 



M 



2^^^-^ 5Z a v<r' C kncr c k'na' 
fi=x,y,z n=l kk'crcr' 

M 



5Z C kn<7 C k'na 



(13) 



n=l cr=T,X kk' 



Here, S is the impurity spin; er M are Pauli matrices; 
J x , J y , and J z are independent spin-exchange couplings; 
Jo/4 represents potential scattering at the impurity site; 
and if is a local applied magnetic field. The Hamilto- 
nian of Eq. (||) corresponds in this notation to M = 2, 
Jx = Jy — 2t) Jz = Jq = 0, and pb9jH = eU [note 
that a labels the channel index in Eq. (^|), whereas in 
Eq. (^) it labels the spin index]. The inclusion of po- 
tential scattering in Eq. (|l^) is mainly intended to allow 
for comparison with the Coqblin-Schrieffer Hamiltonian 
and the self-consistent ladder approximation of Maekawa 
et alB As we shall see, though, its inclusion will provide 
considerable insight into some of the underlying features 
of the NCA. 

We proceed hv introducing Abrikosov's slave-fermion 
representation,c3 which generalizes in__the Anderson 
model to the slave-boson representationO In this repre- 
sentation, one assigns a pseudo fermion to each impurity 
spin state according to 



ft 



I T) , fl 



(14) 



The physical subspace corresponds to the constraint 
Nf = J2 a fafa — 1) w bich represents the fact that the 
impurity spin has only two possible states. This con- 
straint distinguishes the pseudo fermions from ordinary 
fermions. 



4 



A convenient way to calculate physical quantities is to 
work within a grand canonical ensemble with respect to 
the number of pseudo fermions. These are assigned a 
chemical potential A, modifying the Hamiltonian from 7i 
to H — XNf. The projection onto the physical subspace 
is carried out by taking the limit A — ► — oo. Specifically, 
physical averages are given by 

(OW = ^- lim e-^{ONf)x, (15) 

where 

Z imp = lim e-P x {Nf} x (16) 

A — > — 00 

is the "impurity contribution" to the partition function. 
Here subscripts A denote averages within the grand- 
canonical ensemble. In practice, one can drop the Nf 
operator in Eq. (plf ) for those physical operators O that 
give zero when acting on the Nf = subspace, which 
greatly simplifies the calculations. 

A central quantity in the calculation of physi- 
cal observables is the pseudo-fcrmion Green function, 
G^\z), which enters in its projected form G a (z) = 
lim^-oo G„ (z — A). The latter will be refected to 
hereafter as the pseudo-fermion Green function.LJ The 
pseudo-fermion Green function takes the standard form 
G a (z) = [z — e (J — S (T (z)] _1 , where S CT (z) is the projected 
self-energy, and e a = —a^HBgjH is the bare energy of 
the impurity spin (we use u =],[ and er — ± interchange- 
ably to label the spin z component). The Kondo-NCA 
(KNCA) is defined by a particular set of self-energy di- 
agrams, shown in Fig. ^[ Although additional diagrams 
without any crossings do exist for the Kondo-spin Hamil- 
tonian (specifically the particle-particle ladder, see Ap- 
pendix Q), we refer to this approximate scheme as the 
NCA, because it has the same general structure as that 
of the conventional NCA for the Anderson Hamiltonian. 

Working within the slave-fermion representation, we 
write the interaction term of Eq. ( p"3| ) in the form 

M 

YY Y ^(a/3)(a'/3') C nfea C «fe'a'/^'//3) ( 17 ) 
n=l kk' afia'fi' 

where 

V( a /3)( a >f3>) = Y a a a ' a p'f3- (I 8 ) 

fi—0,x,y.z 

Here cr° denotes the 2x2 unity matrix. Viewing the 
^(q/3)(q'/3') coefficients as the elements of a 4 x 4 matrix, 
this choice of labels allows for a convenient summation 
of the ladder diagram of Fig. ||. Details of the derivation 
are provided in Appendix |A|. For a zero magnetic field- 
when both spin orientations are equivalent, one obtainso 

S(e + iS) = y Yl /(- £ fe) A ( £ - c * + ( 19 ) 

k 




k B T/D 

FIG. 3. KNCA results for the zero-field magnetic suscep- 
tibility of the anisotropic two-channel Kondo impurity with 
J z — and different poJ±. Here and throughout the paper 
we use a semi-circular conduction-electron density of states 
with half-width D. At low temperatures, x(^) diverges loga- 
rithmically with decreasing temperature, in accordance with 
exact results. The robustness and universality of the loga- 
rithmic temperature dependence is demonstrated in the in- 
set, where all lines have been collapsed onto a single curve 
by rescaling x(^) an d T with the Kondo temperature, Tk- 
The latter is extracted according to Eq. (^), taking the val- 
ues k B T K /D = 4 x 1(T 4 , 1.15 x 10" 3 , and 3.85 x 10" 3 for 
poJ± = 0.3,0.34, and 0.4, respectively. 

where 

U.(e + iS) = -J2f(tk)G{e + e k +i5). (21) 
fe 

Here /(e) is the Fermi-Dirac distribution function, and 
the coefficients arc defined as 

A = \ (J + Jx + Jy + Jz) , (22) 

At = - (Jo + 2 J ~J X -J y - J z ) (23) 

(i = x,y,z). For Jq = J x = Jy = Jz, corresponding to 
the Coqblin-Schrieffer Hamiltonian, Eqs. (|l9|)-(^l|) co- 
incide with the-self-consistent ladder approximation of 
Maekawa et al£3 

As stated above, Eqs. ([l9|)-(|2l]) are for a zero mag- 
netic field. However, to study the smearing of the charge 
step for the quantum box, it is necessary to include also 
a local magnetic field. The addition of a finite magnetic 
field clearly distinguishes the-present NCA formulati©*, 
from those of the Andersonc.3 and Coqblin-SchriefferEa 
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Hamiltonians. While in the latter two cases one simply 
replaces the zero-field pseudo-fermion propagators enter- 
ing the NCA equations with spin-dependent ones, the 
structure of Eqs. (fl9|)-(|2~i|) is substantially modified once 
H is nonzero. The generalization of Eqs. (19)-(|2l|) to a 
nonzero field is specified in Appendix |a|, Eqs. (|A4| ), (|Aq) , 
and @. 

Although somewhat cumbersome, the magnetic-field 
equations are essential to our discussion. They are solved 
numerically by iterations similar to the zero-field ones. 
From their solution one can compute the field-dependent 
magnetization M(H,T), from which \(H,T) follows by 
direct numerical differentiation with respect to H. In- 
terestingly, we found it advantageous to use this finite- 
field differentiation procedure also when computing the 
zero-field susceptibility. This somewhat surprising result 
stems from the fact that the standard magnetization bub- 
ble acquires a nontrivial vertex correction that is absent 
for both the Anderson and Coqblin-Schricffer Hamiltoni- 
ans, and which contributes an important log-T diverging 
term to the low-temperature susceptibility of the two- 
channel case (see Appendix |A], Figs. |l^ and [ll]). Defer- 
ring all details of the numerical procedure to Appendix [A|, 
below we summarize the main features of the KNCA ap- 
proach. All results presented below were obtained for a 
semi-circular conduction-electron density of states with 
half-width D: p(e) = po\/l — (e/D) 2 , where p is the 
conduction-electron density of states at the Fermi level. 

The KNCA approach presented above has several no- 
table successes. Taking the zerp»temperature limit and 
performing a Muller-HartmannE^I type of analysis, one 
can show that Eqs. (|l^)-(|2l|) reproduce the exact criti- 
cal exponents of the multi-channel (M > 1) Kondo ef- 
fect for antiferromagnetic couplings, similar to the NCA 
treatment of the multi-channel Anderson model.L3 Focus- 
ing on antiferromagnetic couplings, the low-temperature 
susceptibility for the M = 2, two-channel case diverges 
logarithmically with decreasing temperature, in accor- 
dance with exact results.113 For the anisotropic limit of 
Eq. (^) , this is demonstrated in Fig. 0, for different values 
of PqJ± (J± being the transverse spin-exchange coupling: 
J± = J x — J y ). Most significantly, there is quantitative 
agreement with the Bethe ansatz for the field and tem- 
perature dependence of the susceptibility in the isotropic 
limit £3 as shown in Fig. |]. Here and throughout the pa- 
per we extract the two-channel Kondo temperature from 
the Bethe ansatz expression for the slope of— the log-T 
diverging term in the zero-field susceptibility Ej 



X(T) 



{^b9j) 



MTk/T). 



(24) 



20k B T K 

Hence the KNCA correctly describes the non-Fermi- 
liquid fixed point of the multi-channel Kondo effect both 
qualitatively and quantitatively. Finally, Eqs. (p^|)-(|2l|) 
correctly give rise to a quantum phase boundary between 
"antiferromagnetic" (J z > — | | ) and "ferromagnetic" 
(J z < — \J±\) couplings, each of-which is characterized 
by different threshold exponentsxSEll 
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FIG. 4. Temperature and field dependence of the magnetic 
susceptibility within the KNCA: full lines are for the isotropic 
two-channel Kondo model with Jo = and poJ = 0.2667 
{UbTk/D = 2x 10 -3 ); dashed lines are for the anisotropic 
two-channel Kondo model with Jo = J z — and poJ± = 0.4 
(k B T K /D = 3.85 x 10~ 3 ). The applied magnetic field takes 
the values of pBgjH/ksTK = 0,0.1,0.3, I, and 10, according 
to the labels marking each set of curves. For each field, the 
associated temperature T = psgjH /fcg is indicated by an 
arrow. Inset: Comparison between the Bethe ansatz (dotted 
lines, taken from Ref. ^) and the KNCA (full lines), for the 
isotropic two-channel Kondo model. Here poJ = 0.2667 for 
the KNCA curves. The values of psgjH /ksTK used and 
their line assignments are identical to those in the main figure. 



At the same time, Eqs. (|l9|)-([2l|) have several clear 
shortcomings: (i) A spurious Kondo-type effect (i.e., a 
low-tempcrature resonance with an underlying nonper- 
turbative energy scale) is found for ferromagnetic cou- 
plings, albeit with different characteristic exponents than 
in the antiferromagnetic case. For Jo > 0, this spurious 
Kondo-type effect extends also to pure potential scatter- 
ing; (ii) Starting from a particle-hole symmetric Hamil- 
tonian, the KNCA breaks particle-hole symmetry, which 
is manifest, for example, in the asymmetric energy de- 
pendence of the conduction-electron T-matrix [See Ap- 
pendix [A| Fig. (iii) Inaccuracies are found in the 
exponential dependence of Kondo temperature on the 
inverse coupling constants. Specifically, the exponential 
dependence of Tk on 1/ J M can be estimated from the lo- 
cation of the lowest pole in the first iteration for A(e) at 
zero temperature. Subsequent iterations only modify the 
pre-exponential factor. For antiferromagnetic couplings 
with Jj > one thus obtains 



T, 



(NCA) 



K 



D exp 



(25) 



which reduces in the isotropic case [when Aq = Sp^J /4, 



see Eq. (E2)] to T 1 



Dexp| 

conduction-electron bandwidth. 



(NCA) 
K 



-4/3/90 J]. Here D is the 
This should to be con- 



G 



trasted with the well-known result TW-r-' Dexp[— 1/poJ], 
valid for any number of channels Af.E3 Similarly, for the 
anisotropic Hamiltonian of Eq. (||) one obtains x^ CA ^ ~ 
D exp [—2 /po J±] , compared to Tk ~ -D cx p[ — 7r /2po J±] 
of Eq. (E2I). Only in the Coqblin-Schrieffer limit, when 
a nonzero potential scattering Jo = J is included, does 
one recover the correct exponential dependence of Tk on 
1/J for isotropic spin-exchange couplings. 

As discussed at length in Appendix H, we are able to 
pinpoint these shortcomings of the NCA as due to the 
omission of particle-particle diagrams. While particle- 
particle diagrams are of higher order in 1 /N in a formal 
large- N expansion as per the Anderson Hamiltonian,Efl 
they must be treated on equal footing in the spin-^ 
Kondo model. The omission of particle-particle diagrams 
leads to the spontaneous generation of potential scatter- 
ing, drives the system to a ferromagnetic Kondo effect, 
and spoils the exponential dependence of the Kondo tem- 
perature on the inverse coupling constants. It is our belief 
that their omission is also responsible for the well-known 
NCA failure to describe the Fermi-liquid ground state 
of the single-channel Kondo effect. We expand on these 
ideas in Appendix ^|, where an extension of the KNCA 
scheme is proposed based on a summation of all parquet 
diagrams. It is shown that even the crudest treatment 
of the extended scheme gives encouraging results for the 
amendment of the above NCA flaws. 

It is important to note, though, that since Tk is expo- 
nentially sensitive to the microscopic parameters of the 
system, it is best extracted directly from the experiment. 
Thus, viewing Tk as an experimental parameter, one can 
use the KNCA to reliably study the delicate interplay 
between the temperature T, the local magnetic field H , 
and the Kondo temperature Tk in the non-Fermi-liquid 
regime of the two-channel Kondo effect. 
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FIG. 5. Temperature dependence of the charge step be- 
tween the (Q) = —eN and (Q) — —e(N+ 1) charge plateaus, 
as obtained from the KNCA. Here we parameterize the tran- 
sition between the plateaus by 8N — —{Q)/e — N. The tun- 
neling matrix element is taken to be pot = 0.2, corresponding 
to poJ± = 0.4 in the anisotropic two-channel Kondo repre- 
sentation of Eq. ([13]) . The corresponding Kondo tempeipfcure 
is equal to k B T K /D = 3.85 x 10~ 3 . Following Matveev^il we 
identify the conduction-electron bandwidth D with e 2 /Cb, 
such that ell/D is equal to [/Cs/e. The overall shape of 
the charge step is mostly unchanged for T <C Tk, however 
the slope around U = continues to steepen with decreasing 
temperature (see inset). Eventually, an infinite slope develops 
at U = as T -> 0. 

ceptibility of the anisotropic two-channel Kondo Hamil- 
tonian with J z = 0. 



IV. APPLICATION OF THE KNCA TO THE 
COULOMB BLOCKADE 

In this section, we apply the KNCA to study the 
smearing of the charge step for a quantum box, in the 
case of a narrow point contact. Smearing of the charge 
step was previously studied by various authors using 
different approaches. These include renormalization- 
group techniquesxil perturbation theory in the tunnel- 
ing conduct|a*ce,L3 infinite summations of spbclasses of 
diagrams J£JE3 and Monte Carlo simulations E3 However, 
these approaches generally fail at low temperatures near 
the degeneracy point. Here we exploit the KNCA to re- 
liably study this regime for a single-mode junction. 

As discussed in Sec. [n| the excess charge in the quan- 
tum box is related within the two-channel Kondo map- 
ping to the isospin magnetization [see Eq. (|l0|)], while 
the capacitance corresponds to the isospin susceptibility. 
Our task therefore amounts to computing the field and 
temperature dependence of the magnetization and sus- 



A. Evolution of the charge step 

Generally speaking, there are three energy scales that 
govern the shape of the charge step for the quantum 
box: ksT, IzbTk, and e 2 /Cb- At high temperatures, 
ksT ^> e 2 /CB, the classical picture is recovered. Ther- 
mal fluctuations are sufficiently strong to wash out the 
Coulomb blockade, and a linear Q vs. <f> curve is ob- 
tained with a slope of Cb- The Coulomb-blockade pic- 
ture emerge as ksT is lowered below e 2 /Cs- Specif- 
ically, for ksT <C e 2 /Cb there are extended, nearly 
horizontal charge plateaus, separated by rather narrow 
charge steps. The charge plateaus extend over a length 
of A(eZ7) platoau rj e/Cs, while the charge steps arexpn- 
fined to a width of A(eC/) stcp cx max{/c B T, k B T K }Ej It 
is only in this temperature regime and in the vicinity of 
the charge steps that the mapping onto the Hamiltonian 
of Eq. (H) applies. 

Restricting attention to ksT <C e 2 /Cb, Fig. || depicts 
the evolution of the charge step with decreasing T/Tk- 
For Tk <C T, the charge step is governed by thermal 
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FIG. 6. Temperature dependence of the capacitance line 
shape C(U,T). Here the temperatures and line-type assign- 
ments are the same as in Fig. [s| With decreasing temper- 
ature, C(U,T) develops a sharp peak about the center of 
the charge step (U = 0), the height of which diverges loga- 
rithmically with decreasing T. Left inset: The reduced ca- 
pacitance, C(U,T)/C(0,T), plotted versus x x = eU/k B T. 
With increasing temperature, the curves approach the shape 
of the derivative of the Fermi-Dirac distribution function, 
normalized to unity peak height [i.e., 4e a:i /(l + e X1 ) 2 , 
full circles]. Right inset: C(U,T)/C(0,T) plotted versus 
X2 = eU /ks^/TTK, for the four temperatures below Tk (i.e., 
T/T K = 0.0025, 0.025, 0.075, and 0.25). For small to interme- 
diate X2, all lines approximately collapse onto a single curve. 



fluctuations. This is best seen in the capacitance line 
shape, which approaches the shape of the derivative of 
the Fermi-Dirac distribution function for elevated tem- 
peratures (see Fig. [| left inset). For Tk <C T, the width 
and shape of the charge step are therefore set by the tem- 
perature. By contrast, quantum fluctuations dominate 
for T <C Tk- The overall extent of the step is propor- 
tional to IcbTk, and its shape is mostly the same as at 
T = 0. Note, however, that the slope at U — con- 
tinues to steepen with decreasing temperature (see inset 
of Fig. ||), approaching an infinite slope as T — > 0.E1I 
This lack of saturation in the capacitance at U = cor- 
responds to the divergence of the linear susceptibility 
within the two-channel Kondo model, which is a hall- 
mark of the non-Fermi-liquid ground state that develops 
in this case. 

The evolution of the charge step is best resolved in the 
capacitance line shape C(U,T), shown in Fig. [| When 
plotted as a function of UCs/e (corresponding to eU / D 
with D = e 2 /Cb), the capacitance develops a sharp peak 
about the center of the charge step (U = 0), the height 
of which diverges logarithmically with decreasing T. To 
maintain a total integrated weight of one unit charge, 
there is a redistribution of weight from the shoulders of 
the peak to its center as T is decreased. 



The crossover from thermal fluctuations to quantum 
fluctuations can be quantified by considering the reduced 
capacitance, C(U, T)/C(0, T). To this end, we have plot- 
ted the reduced capacitance of Fig. |(| once versus x\ — 
eU/ksT (left inset), and once verses x 2 = eU /IzbVTkT 
(right inset). Clearly, two different qualitative behaviors 
are found, depending on whether T 3> Tk or T 1 C Tk- 
For T Tk, the reduced capacitance approaches the 
curve 4e Xl /(l + e Xl ) 2 , which is just the derivative of the 
Fermi-Dirac distribution function, normalized to unity 
peak height. This is to be expected of thermal fluctua- 
tions, when the occupations of the two available charge 
configurations in the quantum box follow the Fermi-Dirac 
distribution function. A different qualitative behavior is 
seen for T <C Tk- Here there is a strong fanning out 
of the reduced capacitance when plotted versus Xi, but 
an apparent scaling when plotted versus xi- As argued 
below, this approximate scaling with X2 is an important 
characteristic of the non- Fermi-liquid regime of the two- 
channel Kondo effect. 



B. Approximate scaling with U/VT 

To understand the apparent scaling of U with a/T for 
T <C Tk, we go back to the two-channel Kondo represen- 
tation of the Coulomb blockade. As is well known, the 
intermediate-coupling non- Fermi- liquid fixed point of the 
two-channel Kondo model is unstable against an applied 
magnetic field, which drivep-the system to a phase-shifted 
Fermi-liquid ground stateE3 For a weak magnetic field, 
fiBgj\H | <^ Tk, this crossover from non- Fermi-liquid to 
Fermi-liquid behavior is associated with an energy scale 



kBT x = 



{HB9jH)< 
kBT K 



(26) 



that depends quadratically on H. Loosely speaking, for 
Tj < T <C Tk a magnetic field has only a minor ef- 
fect on the low-temperature thermodynamics, while for 
T < T x < Tk the effect of a temperature is small. In 
particular, the susceptibility \{H,T) is weakly field de- 
pendent for T x <C T, and is logarithmically field depen- 
dent for T <§; T x . The connection with the reduced- 
capacitance plot of Fig. ^ is made by recalling that 
eU and C(U, T)/C(0, T) in the Coulomb blockade map 
onto fisgjH and x(H,T)/x(0,T), respectively, in the 
two-channel Kondo representation. Thus, when plotted 
versus X2 — ±\/T x /T, the reduced capacitance shows 
a characteristic crossover from approximately one for 
| X2 1 < 1, to logarithmic in | ar 2 1 for I52I > 1- 

Although the right inset of Fig. || seems to suggest an 
exact scaling form for the reduced capacitance as a func- 
tion of X2 for T <^ Tk, we emphasize that this scaling is 
only approximate. Namely, writing the reduced capaci- 
tance in the form 



F(x 2 ,T) = C(U,T)/C(0,T) 



(27) 
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FIG. 7. Comparison of the proposed scaling function 
F(x2,T) at the Emery-Kivelson point with the one obtained 
at the Coulomb blockade using the KNCA (curves taken from 
the right inset of Fig. ^). Here symbols correspond to the 
Emery-Kivelson point, while lines mark the curves of the 
KNCA. Note that T is varied over two decades of temper- 
ature, from T/T K = 0.25 down to T/T K = 0.0025. 



with X2 — eU /ks^TxT, F does not reduce in the limit 
T -C Tk to a function of X2 alone. Nevertheless, the 
residual T dependence is logarithmic in nature (see be- 
low), and therefore relatively weak. 

So far our discussion has been limited to the KNCA, 
and to the J z — limit of the anisotropic two-channel 
Kondo model. However, we expect the approximate scal- 
ing function of Eq. (E7h to be generic to all antiferro- 
magnetic couplings of the anisotropic two-channel Kondo 
Hamiltonian. To support this conjecture and to gain an- 
alytical insight into the structure of F(x2,T), we con- 
sider below the opposite limit of J z ^> J± 7 for which an 
exact analytic solution of the two-channel Kondo model 
exists for a particular value of J z . Ai_this special point, 
known as the Emery-Kivelson point J£3 one can compute 
F(x2,T) exactly, without resorting to the KNCA. The 
results for this special limit are summarized in Fig. |?] 
and Eqs. below. 



C. Comparison with the Emery-Kivelson point 

As noted by Emery and Kivelson,@ the two-channel 
Kondo model with poJ± <C PoJ z = 1 (Jz = kvf in the 
notations of Ref. p9| ) can be mapped exactly onto a nonin- 
teracting "Majorana" resonant-level model, which is the. 
analog of the Toulouse limit for the single-channel case.EJ 
At this special point, one has an exact analytic solution 
of the anisotropic two-channel Kondo model, which cap- 
tures the essential low-temperature physics of the model. 
In particular, one can use the Emery-Kivelson point to 
extract universal low-temperature behaviors of the two- 



channel Kondo model @ 

Using the exact solution of Emery and Kivelson, we 
computed the local-field and temperature dependence of 
the magnetization at the solvable point. The resulting 
expression reads 

M( g ,T) = ^|^Re{ ^-^ } ) (28) 

where ip(z) is the digamma function,!!! z± are equal to 
1 , $T K 



z± 



and A is given by 



7T 2 T 



(1±A), 



A = Wl- 



10T K 



(29) 



(30) 



Here Tk is the Kondo temperature, defined accord- 
ing to the Bethe ansatz expression for the slope of 
the log-T diverging term in the zero-field susceptibil- 
ity, Eq. With this definition, Tk is related to the 
half-width T of the Majorana resonant level (see Ref. |59| ) 
through Tk — 7rr/20. The field-dependent susceptibility 
%(-ff, T) is readily obtained from Eq. (^8|) by differenti- 
ating M(H,T) with respect to H. 

Focusing on T < T K and substituting C(U, T)/C(0, T) 
and eU for x(H,T)/x(0,T) and ^BgjH, respectively, 
Fig. compares the proposed scaling function F(x2, T) 
at the solvable point with the one obtained within the 
KNCA for the Coulomb blockade. There is good over- 
all agreement between the two cases, even though one 
is dealing with two opposite limits of the anisotropic 
two-channel Kondo model, and with two entirely differ- 
ent methods of calculation (KNCA versus an exact solu- 
tion based on mapping onto a "Majorana" resonant-level 
model). Although the residual T dependence is more 
pronounced at solvable point, it is still quite small con- 
sidering that T is varied in Fig. ^ over two decades of 
temperature. 

The great advantage of the solvable point, though, lies 
in the explicit form of F(x2, T), which allows for an an- 
alytic study of the limit fcgT, e\U\ -C ksTK- Explicitly, 
omitting terms that are either a factor of O ((cU/Tk) 2 ) 
or a factor of O (T/Tk) smaller than the leading-order 
ones, F(x2,T) is expressed as 



F(x 2l T) = 1 



h(T)> 



(31) 



where 



1 x 2 



20 



(i) 



■ r 2 

40 



4, 



(32) 



and 
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JA" = 1 — f(eU) follows the Fermi-Dirac distribution 
function, one has the standard relation 



FIG. 8. The capacitance line shapes of Fig. ^| re-plotted 
as a function of SN = —{Q) /e — N. Within the two-channel 
Kondo representation, this corresponds to plotting the mag- 
netic susceptibility as a function of the magnetization. For 
T < Tk, there are strong deviations from the parabolic form 
of Eq. (|34|), which characterizes pure thermal fluctuations. In- 
set: A plot of the reduced capacitance C(U, T) / C(0, T) versus 
SN. In accordance with predominant thermal fluctuations, 
for T 3> Tk the reduced capacitance approaches the parabola 
4<5AT(1 - SN) (line with full circles). 



2 T I 



"1 5 



log 



7.2T K 
T 



(33) 



Here tp^(z) = dip/dz is the trigamma functionj£3 while 
g(x 2 ) ~ 3ir 2 x?,/80 for \x 2 \ < 1. Clearly, T enters 
F(x 2 ,T) only logarithmically through h(T), which ex- 
plains the rather weak temperature dependence of the 
proposed scaling function. We are therefore led to con- 
clude that a scaling plot of the reduced capacitance ver- 
sus X2 — eU /ks^/TxT provides a sharp diagnostic for 
the two-channel Kondo scenario, directly testing the as- 
sociated non-Fermi-liquid scaling of U with \pT at low 
temperatures. 



D. Capacitance versus average charge 

Finally, a different manifestation of the onset of quan- 
tum fluctuations is provided by the inter-relation between 
the capacitance line shape and the average charge in the 
quantum box. Specifically, in Fig. @ we have re-plotted 
the capacitance line shape of Fig. gas a function of the 
average number of excess electrons in the box, parame- 
terized by SN — —(Q)/e — N. Within the two-channel 
Kondo representation, this corresponds to plotting the 
magnetic susceptibility as a function of the magnetiza- 
tion. In the case of pure thermal fluctuations, when 



C(U,T) 



knT 



SN(1 - SN). 



(34) 



Here we have assumed that k B T,e\U\ -C e 2 /Cs, such 
that all higher energy charge configurations can be ne- 
glected. Thus, up to a prefactor of 1/T, C(U, T) displays 
a universal parabolic dependence on SN. The onset of 
quantum fluctuations is reflected in the breakdown of the 
above parabolic form. 

As seen in Fig. ||, for T Tk the capacitance line 
shape approaches the parabolic form of Eq. (|34|) , in accor- 
dance with predominant thermal fluctuations. By con- 
trast, for T < Tk there are strong deviations from the 
parabolic form of Eq. (j|J) . Specifically, the capacitance 
line shape is much narrower then that of Eq. (|34|), with 
concave shoulders that saturate with decreasing temper- 
ature. This is most clearly seen in the inset of Fig. @ 5 
where the reduced capacitance is plotted versus SN. 



V. DISCUSSION 

In this paper, we have extended the NCA to the 
Kondo-spin Hamiltonian with arbitrary spin-exchange 
and potential-scattering couplings, and applied it to 
the charge fluctuations in a single-electron box at the 
Coulomb blockade. The KNCA correctly describes the 
non-Fermi-liquid fixed point of the multi-channel Kondo 
effect both qualitatively and quantitatively. It repro- 
duces the exact non-Fermi-liquid exponents and loga- 
rithms of the multi-channel Kondo effect, and gives sur- 
prisingly accurate results for the temperature and field 
dependence of the magnetic susceptibility in the isotropic 
two-channel case (see Fig. |[). Hence the KNCA offers a 
reliable approach to study the delicate interplay between 
the temperature T, the local magnetic field H, and the 
Kondo temperature Tk in the non-Fermi-liquid regime 
of the two-channel Kondo effect. 

At the same time, the KNCA reveals several short- 
comings of the NCA that are not apparent from the 
more common applications of this approach to the multi- 
channel Anderson and Coqblin-Schrieffer Hamiltonians: 

1. Kondo temperature — For general antiferromag- 
netic couplings, the NCA class of diagrams is insuffi- 
cient for producing the correct exponential dependence 
of Tk on the inverse coupling constants. In fact, for an 
isotropic spin-exchange interaction J, the correct expo- 
nential dependence of Tk on 1/J is recovered only if a 
matching potential-scattering term of magnitude Jo — J 
is included. Interestingly, this is precisely the case corre- 
sponding to the Anderson and Coqblin-Schrieffer Hamil- 
tonians. That Jo enters the exponential dependence of 
Tk is an artifact of the KNCA, as is the fact that Jo = J 
is required in order to recover the correct exponential 
dependence of Tk on 1/J. 
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2. Particle-hole symmetry - - The NCA inher- 
ently breaks particle-hole symmetry. Starting from a 
particle-hole symmetric Hamiltonian, the KNCA violates 
particle-hole symmetry, which is most clearly manifest in 
the asymmetric energy dependence of the conduction- 
electron T-matrix. We note that a similar asymmetry in 
the shape of the Kondo resonance was also found for the 
multi-channel Anderson model, where it was attributed, 
to the lack of an underlying particle-hole symmetry.E3 
From the KNCA we conclude that such an asymmetry in 
the Kondo resonance is an inherent feature of the NCA, 
independent of whether the underlying Hamiltonian is 
particle- hole symmetric or not. 

3. Ferromagnetic Kondo effect — One regime obvi- 
ously not accessible starting from an Anderson Hamil- 
tonian is that of ferromagnetic spin-exchange couplings. 
For ferromagnetic couplings, the KNCA produces a spu- 
rious Kondo- type effect, which extends to pure potential 
scattering if J > 0. 

As discussed in detail in Appendix [b], all the above 
deficiencies of the KNCA are a consequence of the omis- 
sion of particle-particle diagrams, which must be treated 
on equal footing for the spin-i Kondo model. To this 
end, an extended scheme was devised in Appendix |C], 
based on a complete summation of all parquet diagrams. 
Preliminary studies of the extended scheme seem to sug- 
gest the amendment of the above NCA flaws. We further 
expect it to remedy the well-known NCA failure to de- 
scribe the Fermi-liquid ground state of the single-channel 
Kondo effect, much in the same way as the conserving 
T-matrix approach does for the Anderson model. Al- 
though we anticipate an intimate relation between our 
extended schema, and the conserving T-matrix approach 
of Kroha et al.,c3 we are unable to establish any formal 
linkage between the two approaches. 

Applying the KNCA to the charge fluctuations in a 
single-electron box at the Coulomb blockade, two ap- 
proaches have been proposed to quantify the crossover 
from thermal fluctuations to quantum fluctuations in the 
case of a narrow point contact: (i) By plotting the re- 
duced capacitance C(U,T)/C(0,T) versus appropriate 
scaling combinations of the electrostatic potential U and 
the temperature T; and (ii) By plotting the reduced ca- 
pacitance versus the average number of excess electrons 
in the quantum box, SN. In the latter case, pure thermal 
fluctuations are characterized by a universal parabolic 
form, which breaks down as soon as the temperature is 
lowered down to Tk- In the former case, different scalings 
of U with T characterize the high-temperature (thermal 
fluctuations) and the low-temperature (quantum fluctu- 
ations) regimes. As seen in Fig. |[ for e 2 /Cb T ^> Tk 
the reduced capacitance approaches a universal func- 
tion of x\ — eU/ksT, corresponding to the deriva- 
tive of the Fermi-Dirac distribution function. By con- 
trast, for Tk T there is an approximate scaling with 
X2 = eU/kBy/T~K~T, which differs markedly from both 
thermal fluctuations and ordinary lifetime broadening. 

This unusual scaling of U with y/T for T -C Tk is a dis- 



tinct characteristic of the non- Fermi-liquid regime of the 
two-channel Kondo effect. It directly probes the anoma- 
lous one-half scaling dimension of an applied magnetic 
field near the intermediate-coupling fixed point of the 
two-channel model. Although the reduced capacitance 
does not show exact scaling with x-i , the residual temper- 
ature dependence is logarithmic in nature, and thus suf- 
ficiently weak to render a scaling plot of C(U, T)/C(0, T) 
versus Uf yT a sharp experimental diagnostic for the ob- 
servation of the two-channel Kondo scenario. For exam- 
ple, by plotting the reduced capacitance versus U/VT for 
two temperatures differing by a factor of two or more, one 
should be able to distinguish whether the relevant scaling 
variable is U/y/T or not. The sole requirement is that 
both temperatures be lower than Tk, which itself may 
be extracted from the shape of the scaling curve. More- 
over, the proposed scaling relation should equally apply 
to the case of nearly perfect transmission, which likewise 
corresponds to the two-channel Kondo effect .E3 

A few words are in order at this point about the ex- 
perimental feasibility of reaching the low-temperature, 
non- Fermi- liquid regime of the two-channel Kondo effect. 
There are two main limiting factors for observing a fully 
developed two-channel Kondo effect: (i) The charging 
energy, e 2 /Cb, must be sufficiently large in order for a 
measurable Kondo temperature to emerge; (ii) The mean 
level spacing in the quantum box must be sufficiently 
small as not to cut off the two-channel Kondo effect. 
These two conditions set opposite limitations on the size 
of the quantum box, which cannot be too large Bar too 
small. As recently pointed out by Zarand et al.pl it is 
practically impossible to meet both conditions in present 
semiconducting devices, but these conditions can be ful- 
filled in metallic quantum boxes. The main complication 
with metallic devices has to do with the construction of 
atomic-size junctions, which are required for obtaining 
single-mode junctions. The latter are needed since the 
relevant energy scale for the onset of the two-channel 
Kondo effect decreases exponentially with the number of 
tunneling modes .EJ 

Finally, we wish to mention several potential appli- 
cations of the KNCA, beyond the present study of the 
charge fluctuations in a quantum box. The most natural, 
perhaps, is the application of the KNCA to the nonequi- 
librium scattering off nonmagnetic two-level tunneling 
systems. Previous efforts in this direction have fo- 
cused on the two-channel Anderson modelJ13 which has 
no direct microscopic justification for TLS. The KNCA 
allows one to study the actual anisotropic Kondo model 
appropriate for this case. It would be particularly inter- 
esting to see whether there are any quantitative changes 
to the scaling curves for the differential conductance 
with eV/ksT upon going from the two-channel Anderson 
model to the anisotropic two-channel Kondo model. 

Another interesting possibility is the extension of the 
present work to the calculation of the nonequilibrium 
current through a single-electron transistor near the de- 
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generacy point. In this case, the single-electron box is 
connected to two separate leads, in between which a 
voltage bias is applied. The linear response of such a 
system was cptnputed for strong tunneling by Furusaki 
and Matveevj£-3 under the stringent assumption that each 
lead is coupled to different conduction-electron modes 
within the bps (i.e., no coherent propagation between 
the two leadsEj). This assumption, which corresponds to 
a four-channel Kondo Hamiltonian, has far-reaching im- 
plications for the zero-temperature conductance, which 
vanishes at the degeneracy point for any ^asymmetry in 
the couplings to the left and right leads .E-3 In the gen- 
eral case we expect, though, some overlap between the 
modes coupled to the left and right leads, which should 
restore the two-channel Kondo picture for tunneling at 
sufficiently low temperatures. Whether this results in a 
square-root voltage and temperature dependence of the 
differential conductance as in scattering off TLS remains 
to be seen. 

Lastly, one can exploit the KNCA to study the STM 
spectra around an isolated Kondo impurity placed on 
a e?-wave superconductor. Such a setting may be re- 
alized in recent STM-, experiments on Zn impurities 
in Bi 2 Sr2CaCu 2 8+( 5.n While large- N mean-field the- 
ories for the STM spectra around an isolated magnetic 
adatota .are expected to capture the main qualitative fea- 
turesj23'E3 these approaches are subject to various limi- 
tations, and are quantitatively inaccurate. The KNCA, 
while not free of shortcomings, should provide a more 
accurate framework for describing the STM spectra. 
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APPENDIX A: DETAILS OF THE KNCA 

1. Derivation of the magnetic NCA equations 

In the presence of a nonzero magnetic field, spin up 
and spin down are no longer equivalent. This complicates 
the structure of the KNCA equations, which involve in 
addition to the coefficients of Eqs. (p2|)-(|2^) also the 
following combinations: 



A 



(Jo- 
(Jo 



J,)/4, 
J,)/4, 



r_L 

A± 



(Jx 



(Al) 
(A2) 



Adopting the 4x4 matrix convention of Eqs. Ql'i 
(If), Fig. I shows the building block of the NCA 
self-energy diagrams of Fig. 0, which we denote by 



k net' 
ioo 



(ccPXtf. P' ) 
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FIG. 9. Diagrammatic representation of P^ x \iv), which 
serves as the building block of the NCA self-energy diagram 
of Fig. |^. p( A ) (iu) consists of a particle-hole bubble and a bare 
vertex attached at one end. Here a,/3, a', and /3' are spin la- 
bels, n is a conduction-electron channel index, and there is 
an implicit summation over the conduction-electron momen- 
tum index k. Dashed lines represent fully dressed / Green 
functions. 

0){a',p')- It consists of a particle-hole bubble 
and a bare vertex attached at one end. Upon projection, 
P n (z) = liniA^-co Pn, X \z — A) takes the form 



Pn(z) 



where 



r n T (z) 

AoIFi^) A ± n T (z) 

o A ± n x (z) A n T (z) 

r ± n T (z) o o 



k 



o 
o 

r n x (z) 



(A3) 
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is the projected bubble. Here we have written the matrix 
P n (z) with respect to the basis (t,T) <-> 1, (Jot) *-* 2, 
(t>l) 3, and (J.,!) <-> 4, and made use of the fact 
the pseudo-fermion Green function is diagonal in the spin 
index. Note that P does not depend on the channel index 
n, which is omitted hereafter. 

Using the above matrix convention, the ladder diagram 
of Fig. |2| is conveniently reduced to a geometric series in 
the matrix P. The resulting expression for the projected 
pseudo-fermion self-energy readso 



1 



l-P(z-e k 



■V 



(aer)(acr) 

(A5) 



where M is the number of conduction-electron channels, 
and V is the interaction matrix of Eq. (fl8|) . The final task 
in deriving t he K NCA equations is to invert the 1 — P ma- 
trix in Eq. ( |A5| ). For a zero magnetic field, when Hx(z) 
and Ilx(z) are equal, one readily obtains Eqs. (|l9|)-(pil). 
For a nonzero magnetic field, when Il-f (z) and II; (z) dif- 
fer, a more complicated expression is obtained. Explic- 
itly, E CT (z) takes the form 



(A6) 
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where 
A CT (z) 



-2r + 2(rg-ri)n_ CT (z) 



[l - r n T (z)][i - ron^z)] - rin T (z)n x (z) 

-2A + 2(A2-Ai)n_ ff (z) 



[i - A n T (z)][i - A n x (z)] - A2 n T («)n ; («) ' 



(A7) 



Equations (A4), (|A6|), and (A7) constitute the full set 
of KNCA equations for a general nonzero magnetic 
field. For a zero magnetic field, they properly reduce 
to Eqs. ©-©. 



2. Negative- frequency spectral functions 

Although all relevant physical information is contained 
in principle in the ordinary pseudo-fermion spectral func- 
tions, numerical considerations dictate the introduction 
of yet another set of functions. To see this we note that 
the evaluation of any physical observable necessarily in- 
volves the integral of a pseudo-fermion spectral function 
multiplied by a Boltzmann factor that diverges for large 
negative energies. Hence the dominant contribution to 
such integrals comes from a regime where the pseudo- 
fermion spectral functions are exponentially small. Since 
it is numerically impossible to maintain such an expo- 
nential accuracy when solving the ordinary set of KNCA 
equations, an alternative approach is required. This is 
achieved by absorbing the above-mentioned Boltzmann 
factor into the definition of the negative-frequency spec- 
tral functions, 



a a (e) = --e^°- e hm{G a (e + i5)}, 



(A8) 



and deriving an alternative set of equations directly for 
the functions a CT (e). In the above equation, /? = l/ksT 
is the inverse temperature, and eo is an arbitrary nu- 
merical parameter, chosen to avoid exponentially large 
numbers xB The main physical quantities of interest are 
conveniently expressed in terms of the negative- frequency 
spectral functions. Specifically, the impurity contribu- 
tion to the partition function, Eq. (Uq), is given by 



Z — p~P e " 



M e ) +ai(e)]de, 



(A9) 



while the impurity magnetization is expressed as 



1 f°° fat(e) - ai fell de 

M(H, T) = -^ flg /-°° T ■ (A10) 

2 J-oo M e ) + a l(e)] & 

To derive an alternative set of equations for the 
negative- frequency spectral functions, it is necessary to 
take the imaginary parts of Eqs. (A4), (A6) and (A7), 
and explicitly multiply them by the Boltzmann factor 
e /3(eo-«)^ Upon doing so, one arrives at the following set 
of KNCA equations for the negative- frequency functions: 
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FIG. 10. NCA diagrams for the zero-field susceptibility of 
a Kondo impurity, as derived from the Kubo formula. Here 
filled circles represent dressed vertices (see Fig. while 
dashed lines stand for fully dressed / Green functions. In the 
Coqblin-Schrieffer limit, all higher order terms vanish and 
only the bare bubble is left. In contrast, all terms contribute 
to the susceptibility of a generic Kondo Hamiltonian. 
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Here v?o-(e) and Vv(e) are t w0 rather complicated expres- 
sions given by 



2C 2 (r ,r ± ,n_ 



rin T nj 



(i-ron T )(i-r nj 

+ (the same with To, T 'j_ and Ao, Aj_ interchanged) 
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2C(r ,r ± ,n_ CT )C(r ,r ± ,n, 
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(i-r n T )(i-r n i )-rin T n i 

+ (the same with Tq, T± and Ao, Aj_ interchanged) 
where 



C(x, y,z) = x + (y 2 - x 2 )Re{z}. 



(A15) 



In Eqs. (|AT^)-(|A14|), we have omitted for conciseness the 
arguments of (e + iS) and (e + iS) . 

Note that the negative-frequency equations explic- 
itly depend on G a (z) and H a (z), hence the negative- 
frequency and ordinary KNCA equations must be solved 
simultaneously. This is done by repeated numerical it- 
erations until convergence is reached. As is always the 
case with numerical solutions of NCA-type equations, the 
key to high-precision numerics is in a well-designed grid 
of mesh points that scatter mesh points more densely 
near the threshold and peaks of the relevant functions. 
To this end, we have used a combination of linear and 
logarithmic grids. As a critical test for the precision of 
our numerical code, we have checked in all our runs that 
the pseudo-fermion spectral functions fulfilled the spec- 
tral sum rule to within one part in one thousand (i.e., 
0.1%). 



13 



CD 



o 



300 



200 - 



100 



Anisotropic Kondo: Full 

♦ Anisotropic Kondo: Bubble 
Isotropic Kondo: Full 

i Isotropic Kondo: Bubble 
Coqblin-Schrieffer: Full 

• Coqblin-Schrieffer: Bubble 




k B T/D 



FIG. 11. Comparison between the simple magnetization 
bubble and the full zero-field susceptibility, obtained by dif- 
ferentiating M(H,T) with respect to H. All curves refer to 
the two-channel case, with pnAn held fixed at 0.2. While the 
two procedures give identical results in the Coqblin-Schrieffer 
limit (we use prjJo = poJ = 0.2), substantial deviations are 
found for the isotropic and anisotropic Kondo models (we use 
PoJ = 0.2667 for the isotropic Kondo model, and J z = and 
poJ± ~ 0.4 for the anisotropic Kondo Model). In the latter 
two cases, the effect of the magnetic vertex correction is to 
significantly reduce the slope of the log-T diverging term in 
X(T). 



3. Magnetic vertex correction 
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FIG. 12. Imaginary part of the conduction-electron 
T-matrix, as obtained from the KNCA. Here, the solid line 
corresponds to the anisotropic two-channel Kondo model 
with J z = and p J± = 0.4 (k B T K /D = 3.85 x 10" 3 ); 
the long-dashed line is for the isotropic two-channel Kondo 
model with p J = 0.2667 {k B T K / D = 2 x 10~ 3 ); and the 
short-dashed line is for the two-channel Coqblin-Schrieffer 
model with p J = p J = 0.2 (k B T K /D = 1.35 x 10~ 3 ). The 
temperature in all curves is equal to T/Tk = 0.01. In ac- 
cordance with the analytic treatment of the NCA equations 
at zero temperature] 3 -!! the T-matrix has a cusp at zero en- 
ergy, with a strong particle-hole asymmetry. All curves fall 
within 1.5% from the predicted zero-temperature NCA value 
of 2np Im{T kk , (0 - id)} = tt 2 /8 (marked by an arrow in the 
inset). 



As mentioned in the main text, we have calculated the 
magnetic susceptibility x(H,T) b y num erically differen- 
tiating the magnetization of Eq. (A1C) with respect to 
H. An interesting observation has to do with the cal- 
culation of the zero- field susceptibility %(T). For the 
Anderson and Coqblin-Schi'icffer models, one tradition- 
ally computes the zero-field susceptibility from the Kubo 
formula- which reduces in the NCA to a simple bubble di- 
agram.EJ For a Kondo model, however, there is a nontriv- 
ial magnetic vertex correction that cannot be neglected. 
Fig. [10| depicts the NCA susceptibility diagrams for a 
generic Kondo model. In addition to the simple bub- 
ble, there is a magnetic vertex correction in the form of 
a ladder diagram. Each rung of the ladder consists of 
two dressed vertices (displayed in Fig. |ll|), connected by 
a conduction-electron bubble. Since the rungs are fre- 
quency dependent, a direct evaluation of the entire series 
requires the solution of an integral equation, which is far 
more complicated than numerical differentiation of the 
magnetization curve. 

While in the Coqblin-Schrieffer limit, Jq = J x = 
Jy = Jz, it Possible to show that only the simple bub- 
ble survives ,E2l for a generic Kondo model there are non- 
negligible contributions coming from the magnetic ver- 



For the two-channel case, this is demon- 
11. where the results of the bare bubble 



tex correction, 
strated in Fig. 
are compared with the full zero-field susceptibility ob- 
tained by differentiating M(H,T). As required, the two 
procedures are indistinguishable in the Coqblin-Schrieffer 
limit, establishing the numerical accuracy of our differen- 
tiation routine. A different picture is recovered, however, 
for the isotropic and anisotropic Kondo models. In both 
cases the slope of the log-T diverging term in x(T) is sig- 
nificantly reduced upon going from the simple bubble di- 
agram to the full susceptibility. This effect is particularly 
pronounced for the anisotropic Hamiltonian of Eq. 
where nearly a three-fold reduction is found in the slope 
of the log-T component of x(T). Thus, the magnetic 
vertex correction contains essential contributions to the 
low-temperature susceptibility of the two-channel Kondo 
model. 



4. Conduction-electron T-matrix 

Another quantity of interest in the theory of dilute 
magnetic alloys is the conduction-electron T-matrix, de- 
fined by the expansion 
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r (z) = G<®{z)5 ktk r + G(°)(z)%, ff (2)4°» (*). 
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Here Gkk',a(z) is the Fourier representation of the dressed 
conduction-electron Green function 



Gkk'Ar) = -(T T c fe „ CT (T)cL„ CT (0)), 



(A17) 



and G k °^(z) = l/(z— e&) is the bare conduction-electron 
Green function. For a nondegenerate Anderson impu- 
rity, Tkk',a(z) is equal to V 2 Gi mp . a (z), where V is the 
hybridization matrix element (assumed here to be k in- 
dependent), and Gi mPt o-(z) is the dressed impurity Green 
function. Hence the calculation of the T-matrix is re- 
duced to that of the dressed impurity Green function. 
For a Kondo impurity there is no analogous relation, and 
Tkk',a(z) needs to be constructed directly from Eq. (Alt). 



Within the slave-fermion representation, the T-matrix 
is related to the unprojected conduction-electron self- 
energy through 



T 



kk'a 



00 = 
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Using the NCA class of diagrams for the conduction- 
electron self-energy, and in the absence of an applied 
magnetic field, one obtains 



T k k' (z) 



-0e o 



[d(e)G(z + e) - o(e)A(e - z)] de, 

(A19) 



where G(z) is the pseudo-fermion Green function, q(e ) 
is the negative- frequency spectral function of Eq (|A£ ) , 
A(z) is defined in Eq. (|20|), and d(e) is given by Eq. ( |A12 ). 
Note that all spin indices have been dropped in Eq. ( A19| ) 
due to the equivalence of the two spin orientations for 
a zero magnetic field. Also notice that Tkk 1 (z) has no 
explicit dependence on k and k\ as is characteristic of 
purely local scattering. 

Figure [jj] shows the imaginary part of the KNCA 
conduction-electron T-matrix, for each of the Coqblin- 
Schrieffer, isotropic Kondo, and anisotropic Kondo mod- 
els. The coupling constants were adjusted in each model 
as to give poA — 0.2, while the temperature T/Tk — 
0.01 is sufficiently low in order for fully developed Kondo 
peaks to be seen. Two points are noteworthy. First, 
there is a clear asymmetry in the energy dependence of 
Im{Tcfc' (e + iS)}, which for the isotropic and anisotropic 
Kondo models is in violation of the underlying particle- 
hole symmetry. This asymmetry persists from high ener- 
gies, |e| ~ D, down to at low energies, |e| < Tjc- As 
discussed in Appendix it stems from the omission 
of particle-particle bubbles within the NCA. Thus, the 
asymmetric line shape of the Kondo resonance is an in- 
herent feature of the NCA class of diagrams, and is not 
due to the lack of an underlying particle-hole symme- 
try, as previously asserted for the multi-channel Ander- 
son model.Ba 
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FIG. 13. Renormalization of the bare interaction within 
the KNCA. In the diagrams of Fig ^, the bare vertex is 
dressed by a ladder of particle-hole bubbles, each composed 
of one (bare) conduction-electron line and one (fully dressed) 
pseudo-fermion line. 



The second point to notice is that all curves ap- 
proach the predicted zero-teaaperature NCA value of 
27T/9oIm{T/ £ fc/ (0 — iS)} = 7r 2 /8.c3 Since this value is deter- 
mined by the NCA threshold exponents, which are un- 
changed in going from the two-channel Anderson to the 
isotropic and anisotropic Kondo models, then the same 
zero-temperature value applies to all models. 



APPENDIX B: COMPARISON WITH 
POOR-MAN'S SCALING 

Considerable insight into the origin of the NCA set- 
backs can be gained by corap,aring the KNCA with An- 
derson's poor-man 's scaling.O Specifically, the KNCA set 
of diagrams can be viewed as a ladder renormalization of 
the bare interaction, as illustrated in Fig. [ll| The ba- 
sic renormalization of the bare vertex is due to particle- 
hole bubbles. This should be contrasted with poor-man's 
scaling, where two distinct types of renormalizations are 
found: particle-hole and particle-particle bubbles (see 
Fig. ||). Thus, it appears that the KNCA corresponds 
to retaining only the particle-hole processes within poor- 
man's scaling (although additional higher-order terms do 
enter the KNCA through the self-consistent dressing of 
the pseudo-fermion Green functions). 

To substantiate this correspondence between the 
KNCA and poor-man's scaling, we have rewritten the 
well-known renormalization-group (RG) equations for 
the dimensionless coupling constants J M = pqJ^, (// = 
0, x, y, z), by separating the contributions of the two di- 
agrams of Fig. 14, Labeling the contributions of the 
particle-hole and particle-particle diagrams by ph and pp, 
respectively, the corresponding equations read 



and 



dJ l 
~dT 

dJp 
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dJj 
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JjJk + JqJi 
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FIG. 14. Renormalization of the interaction within 
poor-man's scaling. We distinguish between (a) processes 
involving the excitation of a high-energy electron (parti- 
cle-hole bubble) and (b) processes involving the excitation 
of a high-energy hole (particle-particle bubble). The KNCA 
retains only the particle-hole bubble of (a), albeit with 
self-consistent dressing of the pseudo-fermion Green func- 
tions. 



dJp 
dl 



"4 l2i 

p=0,x,y,z 
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where I is equal to ln(Do/D), and i,j, and k represent a 
cyclic permutation of x,y, and z. Here D is the running 
bandwidth, and Do is the bare bandwidth. By combining 
the two sets of equations, listed above, one recovers the 
standard RG equations£il 



dl ~ JjJk > 
dJ 



dl 



= 0. 



(B5) 
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If we artific ially i gno re the particle-particle contribu- 
tions of Eqs. @-(g) and settle with Eqs. @-@) 
for dJ^ /dl, it is easy to verify that 



A = - (j + J X + Jy + J 2 ) , 

A\ = \ (Jo 



2Jj J x Jy Jz 



(i = x, y, z) obey the single differential equation 



dA„ 



dl 



- A 2 



(B7) 
(B8) 

(B9) 



ph 



Equation (B9) has the exact same form as the stan- 
dard RG equation for the dimensioeJess exchange cou- 
pling in the isotropic Kondo modelFJ As is well known, 
its outcome depends on the sign of the bare coupling 
constant. If is initially positive ( "antiferromagnctic" 
coupling) then this coupling flows towards strong cou- 
pling; otherwise ("ferromagnetic" coupling) it flows to 
weak coupling. In the former case, a nonperturbative en- 
ergy scale T M ~ Z?oe x P [ — 1/poAu] emerges, below which 
A, of Eq. (B9) becomes of order unity. If several A^ 



simultaneously flow towards strong coupling, then the 
largest of the T M 's will emerge as the actual low-energy 
scale. 



From this analysis it is clear that Eqs. (B1)-(B2) flow 




FIG. 15. RG flow trajectories corresponding 
to Eqs. (B1)-(B2) with an isotropic spin-exchange coupling 
J. There are four different RG regimes: (a) — 3J < Jo < J, 
where only Ao is relevant; (b) J < Jo < — 3J, where only Ao 
is irrelevant; (c) — 3J, J < Jo, where all the A M 's are relevant; 
and (d) Jo < — 3J, J, where all the A^'s are irrelevant. The 
flow in region (a) and in the J > part of region (c) is towards 
Jo, J = 00 and Jo/J = 1; in region (b) and in the J < part 
of region (c) the flow is towards Jo, — J = 00 and Jo/J = —3; 
the flow along the J = line in region (c) is towards Jo = 00 
and J = 0; and in region (d) the flow is towards Jo = J = 0. 



of Eqs. (p2|)-(p^) is positive. The underlying nonpertur- 
bative energy scale that emerges then is given by 



p(ph) 
K 



D exp(- 



1 



PoA n 
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towards strong coupling whenever at least one of the A^s 



where A max is the largest of the positive A^'s. Note that 
for antiferromagnetic couplings with J± > 0, j4 max is 
always given by Ao- 

An identical picture is recovered within the KNCA 
equations for a zero magnetic field, Eqs. (|l9|)-(^l|). Here 
the flow towards strong coupling is manifest in the divefc. 
gence of A(e + iS) at the threshold energy as T — > 0,E3 
which happens whenever at least one of the A^s is posi- 
tive. If more than one A^ coefficient is positive then only 
that A^/W — AfJI(e + iS)] component of A(e + id) cor- 
responding to-* dmax diverges, while all other components 
remain finite.LJ The nonperturbative energy scale under- 
lying the NCA is then precisely given by Eq. (BIO), estab- 
lishing the intimate connection between the KNCA and 
the restricted version of poor-man's scaling where only 
particle-hole processes are retained. (Obviously, the two 
approaches are not entirely equivalent, as the KNCA also 
includes self-consistent dressing of the pseudo-fermion 
Green functions.) 

We can now exploit this basic relation between the 
KNCA and Eqs. pl])-(p2]) to clarify the setbacks of 
the KNCA. For simplicity, let us focus on the isotropic 
case, J x = J y = J z = J. From Eq. (B2) one read- 
ily sees that, starting from Jo = 0, a nonzero poten- 
tial scattering is immediately generated. This is to be 
contrasted with the full set of RG equations, in which 
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potential scattering is exactly marginal. Thus, start- 
ing from a particle-hole symmetric Hamiltonian, the 
omission of the particle-particle contributions necessar- 
ily breaks particle-hole symmetry, leading, for example, 
to the asymmetric energy dependence of the conduction- 
electron T-matrix within the KNCA (see Fig [l2|). 

More significantly, the interplay between the spin- 
exchange couplings and the ever-increasing potential 
scattering substantially alters the RG flow diagram. 
In Fi g, [lg , the RG flow trajectories corresponding to 
Eqs. (ftlfl-p^) are plotted in the J — Jq plane. While 
the full set of RG equations, Eqs. (p5|)-(p6|), support 
only hori zo nta l trajectories, non of the trajectories of 
Eqs. ( |Bl| ) (B2) arc horizontal. Overall there are four 
different regimes: (a) — 3J < Jo < J, where only Aq is 
relevant; (b) J < Jq < — 3J, where only Aq is irrelevant; 
(c) — 3J, J < Jo, where all the A^s are relevant; and (d) 
Jo < —3 J, J, where all the A,, 's are irrelevant. Thus, 



with the exception of region (d), Eqs. (B1)-(B2) flow to- 
wards some kind of strong coupling, including in the case 
of ferromagnetic coupling. The emergence of a ferromag- 
netic Kondo-type effect within the KNCA is therefore an 
artifact of the omission of particle-particle diagrams. 

Finally, if one naively extends Eqs. (B1)-(B2) to the 
strong-coupling regime, then four different fixed points 
are recovered: (i) Jo, J = oo and Jq/ J — 1 [approached 
from region (a) and from the J > part of region (c)]; 
(ii) Jo, — J = oo and Jq/ J = —3 [approached from region 

(b) and from the J < part of region (c)]; (hi) Jo = oo 
and J = [approached from the J = line in region 

(c) ]; and (iv) Jo = J = [approached from region (d)]. 
Interestingly, from the solution of the KNCA equations 
at T = one also finds four different low-temperature 
behaviors, corresponding to the precise four basins of at- 
traction listed above. 



APPENDIX C: BEYOND THE NCA 

From the discussion in the preceding Appendix it is 
clear that any proper theory of the Kondo effect must 
include particle-particle and particle-hole processes on 
equal footing. Diagrammatically, this necessitates going 
beyond the NCA class of diagrams, which is the objective 
of the present Appendix. 

Figure |l^ displays a formally exact representation of 
the pseudo-fermion self-energy. The first line in Fig. [l6| 
represents all pseudo-fermion self-energy diagrams. The 
full circle marks the full vertex function, by which we 
mean the collection of all connected four-legged diagrams 
with one incoming and one outgoing pseudo-fermion line, 
and one incoming and one outgoing conduction-electron 
line. We classify the vertex diagrams according to two 
classes. The particle-hole irreducible (PHI) vertex func- 
tion is defined as the sum of all vertex diagrams that 
cannot be disconnected into two distinct vertex diagrams 
by cutting one internal conduction-electron line and one 
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FIG. 16. A formally exact representation of the 
pseudo-fermion self-energy in terms of irreducible vertex func- 
tions. The first line represents all pseudo-fermion self-energy 
diagrams. The full circle marks the full vertex function, i.e., 
the collection of all connected four-legged diagrams with one 
incoming and one outgoing pseudo-fermion line, and one in- 
coming and one outgoing conduction-electron line. The full 
vertex function is expressed in turn in the second line via 
the more elementary particle-hole irreducible (PHI), parti- 
cle-particle irreducible (PPI) , and both particle-hole and par- 
ticle-particle irreducible (PHI&PPI) vertices. For precise def- 
initions, see text. 



internal pseudo-fermion line propagating in opposite di- 
rections (third line of Fig. |l6|). The particle-particle 
irreducible (PPI) vertex function is defined by a simi- 
lar collection of diagrams, that cannot be disconnected 
into two distinct vertex diagrams by cutting one internal 
conduction-electron line and one internal pseudo-fermion 
line propagating in the same direction (fourth line of 
Fig. lsj) . Note that each vertex diagram is a member 
of at least one of the two classes PHI and PPI. Some ver- 
tex diagrams, like the bare vertex, are both PHI and PPI. 
We refer hereafter to these joint diagrams as PHI&PPI. 
Hence the full vertex function is equal to the sum of the 
PHI and PPI vertex functions, minus the PHI&PPI dia- 
grams which are counted twice. 

There are infinitely many diagrams that are PHI&PPI, 
the first two of which are depicted on the fifth line of 
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Fig. |l6|. Provided one is able to sum the entire class of 
PHI&PPI diagrams, then Fig. [if] constitutes an exact 
representation of the pseudo-fermion self-energy, as well 
as of the PHI, PPI, and full vertex functions. In prac- 
tice, it is impossible to sum all the PHI&PPI diagrams. 
However, any selection of a subclass of PHI&PPI dia- 
grams that preserves particle-hole symmetry defines an 
approximation scheme that necessarily respects particle- 
hole symmetry. The simplest choice for the PHI&PPI 
vertex function is to settle with the bare vertex, which 
defines the parquet class of diagrams. Below we focus on 
this particular choice. 

Each of the PHI, PPI and full vertex functions de- 
pend on three energy variables, as well as on the spin 
indices of the four external legs. Adopting the nota- 
tion of Eq. Jl7|), we label the spin index of the outgo- 
ing (incoming) conduction-electron line by a {a'), and 
that of the outgoing (incoming) pseudo-fermion line by 
(3' (f3). As for the energy variables, we choose to work 
with the following combinations: (i) The particle-hole en- 
ergy, Zh, corresponding to the outgoing pseudo-fermion 
frequency minus the incoming conduction-electron fre- 
quency; (ii) The particle-particle energy, z p , correspond- 
ing to the incoming pseudo-fermion frequency plus the 
the incoming conduction-electron frequency; and (iii) 
The pseudo-fermion energy, zf, corresponding to the out- 
going pseudo-fermion frequency. For each of the PHI, 
PPI, and full vertex functions, the projected functions 
are defined as 

Fa/3;a'0'(Zh, z pi z f) = 

, lim T aln'0>( Z h- X,Z p -X,Zf- A), (CI) 

where F^},g, is the corresponding unprojected function. 

Approximating the PHI&PPI vertex function by the 
bare vertex, one can show that the projected PHI, PPI, 
and full vertex functions have the form 

T {PHI) (z h z zA - 

J E H ^ Z h^V^ Z f) a a a ' a $'^ ( C2 ) 



fi= 0,x,y,z 
^aP;a'l3'( Zh i Z P> Z f) ~ 



fi=0,x,y,z 

r {Full) ( Zh z zA - 

1 a/3;a'/3'V Z h, Z P> z f ) " 



2 E M z h,Zp,Zf)^ a ,a^ . 



(C3) 



(C4) 



fi—0,x,y,z 



Thus, each of the above vertex functions retains the form 
of the bare interaction, only with renormalized couplings 
that are energy dependent. Furthermore, one has the 
identity J, — TL^ + — J M , where J M are the bare 
coupling constants. The PHI and PPI vertex functions 
are given in turn by the following set of coupled integral 
equations: 



V (z h , z p , z f ) 
Vi(z h , z p , z f ) 

H (z h , z p , z f ) 
Hi(z h ,z p ,Zf) 



Jo 

Ji- 



= ^ 







■ [Ho © Ji\ - [Hi © Jo] 

y t E ijmV~Lj © i7m]j 
l,z 

[P»eJ^ 



ij'tn L' "J 

j,m— x,y,z 

*' E 

/j.— 0,x,y,z 
J t + [Po © Ji] + [Pi © Jo] 

j,m— x,y,z 



(C5) 

(C6) 
(C7) 

(C8) 



where tij m is the antisymmetric tensor, i is equal to x, y, 
or z, and the convolution operators © and © stand for 



1 



•Ju] = 2 J~] f( e k)G(z h + e k )H^(z h , z p + Ai, zf) 
k 



xJ l ,(z h ,Zp + A 2 ,Zf + A 2 ), (C9) 

[P» eJ v ] = z Y, f(- e k) G ( z P - tk)JA z h - a 2 , z p, z f 

k 

xV^z h -A 1 ,z p ,Zf-A 1 ). (CIO) 
Here Ai and A2 are shorthands for 



Ai 

A 2 



z f 
z h 



z f- 



(Cll) 
(C12) 



(Note that the frequencies z p — Zf and Zf — Zh are just the 
outgoing and incoming conduction-electron frequencies, 
respectively.) Omitting the constant term (AlJoN e )/A 
(N e being the total number of conduction electrons), the 
pseudo-fermion self-energy is equal to 



£(*) 



AI 



E 

fe,fc' 



fi=0,x,y,z 



f(e k )f(-e k AG(z + e k ~e kl )x 
Jf,J^( z -^k', z + ^k, z )- (C13) 



Equations (|C^)-(|C1^) represent a complete summa- 
tion of the parquet class of diagrams, which includes 
in particular all diagrams contained within poor-man's 
scaling. While we believe these equations are intimately 
related, to the conserving T-matrix approach of Kroha 
et alt2 for the multi-channel Anderson model, we are 
unable to formally link the tw o ap proaches. We fur- 
ther emphasize that Eqs. (C5)-( |C13|) apply to a general 
anisotropic Kondo model, ferromagnetic or antifcrromag- 
netic, and hence extend beyond the Schricffcr- Wolff limit 
of the Anderson Hamiltonian. 



A full self-consistent solution of Eqs. (C5)-(|C13j) is a 
difficult task, due to the three energy variables enter- 
ing the coupled integral equations for Vp,[z,h, z p, z f) and 
7~tfj,( z h, z p, z f)- Below we resort to a crude approximate 
treatment, aimed at demonstrating the amendment of 
the KNCA flaws. 
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Since the onset the Kondo effect is manifest in the 
development of strong scattering at the Fermi level, we 
set the energies and thus frequencies of the incoming and 
outgoing conduction electrons equal to zero. [We empha- 
size that, upon implementing the projection of Eqs. (CI), 
the conduction-electron frequencies acquire the value of 
the bare conduction-electron energies.] We begin by not- 
ing t hat t he energies Ai and A2 in each of Eqs. (C9) 
and (CIO) correspond to the single-particle energy dif- 
ference between the conduction electrons incoming and 
outgoing one of the vertex function in that equation — 
J v for Ai, and or for A 2 . Thus, since scattering 
in the Kondo problem is strong near the Fermi level, one 
expec ts t he ma in contribution to the k summations in 
Eqs. (C9)-( |C10| ) to come from that range in k where Ai 
and A2 are small. As a first app rox i mation , we there- 
fore decouple the sums in Eqs. (C9)-( C10 ) by setting 
Ai = A 2 = 0. 

The second approximation we adopt is to linearize 
Eqs. (|Cg|)— (|CT0|) by replacing the PHI and PPI vertex 
functions, H.^ and Vn, with the bare vertex, J„. This 
clearly is a crude approximation, motivated primarily 
by technical considerations. It encompasses, however, 
the ladder summation of the KNCA. Using the above 
two approximations in combination with the identity 
Jii = V^+Ti,^ — J^, we arrive at the following set of equa- 
tions for the full vertex functions = J^iz^z^zf): 

jo = Jo+\m+(zf)-n-(z,)] E j a ( C14 ) 

fi—0,x,y,z 

J l = J l + \ [IL+(zf)-n_(z f )] [JoJi+JiJo] 

+ i[n+(z / )+n_(z / )] J2 £ ijmJjJm- (C15) 
j,m— x : y ,z 

Here, as before, i is equal to x, y, or z, while Tl±(z) are 
the projected particle-hole and particle-particle bubbles: 



n±(z) = -J2f(±tk)G(z±e k ). 



(C16) 



Similar to the KNCA, one can estimate the expo- 
nential dependence of the Kondo temperature on 1/J M 
by substituting the bare pseudo-fermion Green func- 
tion into Eq. (C16), and examining the location of 
the zero-temperature poles in the full vertex functions 
J^(zf,Zf,Zf). For a particle- hole symmetric band, when 
n + (z) and n_(z) coincide, J (zf,Zf,Zf) is simply equal 
to Jq. Thus, there is no renormalization of the potential 
scattering at the Fermi level, in accordance with poor- 
man's scaling. In the isotropic case, J x = J y = J z = J, 
the remaining components of the full vertex function are 
also isotropic, and are given by 



J 



1 - JU(z f ) 



(C17) 



the absence of a ferromagnetic Kondo effect. By con- 
trast, J(zf,Zf,Zf) diverges for J > and T = at 
e r~j — _Dexp[— X/poJ], which we identify with — fc^Tx. 
Hence the correct exponential dependence of Tk on 
1/J is properly recovered. For the anisotropic Hamil- 
tonian of Eq. (||), a similar analysis gives ksTx ~ 
D exp[— \/2/poJ], which still deviates from the desired 
form of Tk ~ Dexp[— Tt/2poJ±], but marks a significant 
improvement over the KNCA. 

Although the above treatment is too crude to convinc- 
ingly establish the amendment of the NCA flaws, we take 
it as an encouraging indication for the elimination of the 
spurious ferromagnetic Kondo effect and the correction 
of the exponential dependence of Tk on l/J^ within the 
proposed scheme. 



Setting T — and Zf = e + iS, one sees that there 
are no divergences in J{zf,Zf,Zf) for J < 0, implying 
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